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Abstract 

This paper concerns the asymptotic expansion of the solution of the Dirichlet-Laplace problem in 
a domain with small inclusions. This problem is well understood for the Neumann condition in di¬ 
mension greater or equal than two or Dirichlet condition in dimension greater than two. The case of 
two circular inclusions in a bidimensional domain was considered in 1^. In this paper, we generalize 
the previous result to any shape and relax the assumptions of regularity and support of the data. Our 
approach uses conformal mapping and suitable lifting of Dirichlet conditions. We also analyze config¬ 
urations with several scales for the distance between the inclusions (when the number is larger than 2). 
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1 Introduction 

In many application fields ranging from elecfrical engineering fo flow around obsfacles, one has fo 
consider Poisson equation in a domain presenting holes. The presence of small inclusions or of a 
surface defecl modifies fhe solufion of fhe Laplace equation posed in a reference domain Dq. If the 
characteristic size of the perturbation is small, then one can expect that the solution of the problem 
posed on the perturbed geometry is close to the solution of the reference shape. Asymptotic expansion 
with respect to that small parameter -the characteristic size of the perturbation- can then be performed. 

The case of a single inclusion a;, centered at the origin 0 being either in Qq or in 5Do, has been 
deeply studied, see [141 [1511201 1711^. The techniques rely on the notion of profile, a normalized 

solution of the Laplace equation in the exterior domain obtained by blow-up of the perturbation. It 
is used in a fast variable to describe the local behavior of the solution in the perturbed domain. For 
Neumann boundary conditions in dimension greater or equal than two and Dirichlet boundary condition 
in dimension greater than two, convergence of the asymptotic expansion is obtained thanks to the decay 
of the profile at infinity. 

The case of some inclusions was considered for example in the series of papers of A. Movchan 
and V. Maz’ya lU^lTTll where an asymptotic approximation of Green’s function is built and justified 
in a domain wifh many inclusions. The poinfs where fhe inclusions are shunk are fixed in fhose works. 
In [2, the Neumann case where the distance between the holes tends to zero but remains large with 
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respect to their characteristic size was investigated for two perfectly insulated inclusions: a complete 
multiscale asymptotic expansion of the solution to the Laplace equation is obtained in a three scales 
case. 

Recently, V. Bonnaillie-Noel and M. Dambrine have considered in ||T] the case of two circular 
defects with homogeneous Dirichlet boundary conditions in a bidimensional domain. They distinguish 
the cases where the distance between the object is of order 1 and the case where it is larger than the 
characteristic size of the defects but small with respect to the size of the domain. They have derived 
the complete expansion and built a numerical method to solve the problem. Our aim is to extend their 
result to any geometry for the inclusion as well as to richer geometric configurations. Note that L. 
Chesnel and X. Clayes have proposed an alternative numerical scheme in Q. 

Let us make precise the problem under consideration. Let Oq be a bounded simply connected 
domain in We consider N inclusions of size e: for any i = 1,... , N let ul := + ew* where w* 

is a bounded simply connected open set containing 0. The centers of the inclusions x* are distinct, 
belong to flo and tend to Xq G flo as e —0. We allow to have the same limit position Xq = Xq for 
some i ^ j but we assume that the distance for any e > 0 is much larger than the size of the inclusion 
(see (14.11) in Section |4] for precise definitions). Therefore, for e small enough, we have w* C flo and 
cj® n a;| = 0 for any i / j. The domain depending on e is then : 

N 

Os := ffo \ ( U ■ 

i=l 

We search an asymptotic expansion for the solution Us of the Dirichlet-Laplace problem in 0^: 

f = f in Os, 

\ Us = 0 on c?Os, 

where f G H~^+^(Oo) with ^ > 0. Since the H^-capacity of w® tends to 0 as e tends to 0 (see e.g. ifT^ 
for all details on the capacity), it is already clear that converges strongly in Hg(Oo) to uq which is 
the unique solution of the Dirichlet-Laplace problem in Oq: 


J —Auo = f in Oq, 

\ uo = 0 on OOq. 


( 1 . 2 ) 


In this convergence, we have extended by zero inside the inclusions. At order 0, we thus have 
Us = uo + Xg where ||rg ||hi(Oq) —)• 0 as e —)• 0. Let us notice that the remainder satisfies 


—Ar® = 0 in Os, 

< Xg = 0 on OOo, (1-3) 

^Xg = —uo on dul, 's/I < i < N. 

The purpose of fhis paper is fo gef fhe lower order ferms in fhe expansion of u^. We follow fhe leading 
ideas of lUl and generalize fheir resulfs fo any geomefry of fhe defecfs while simplifying fhe proofs. 
We also relax fhe assumpfion on fhe regularify of f and do nof more assume fhaf f is supported far 
away fhe inclusions. The drawback is fhaf fhe correctors builf in fhis work are less explicif buf remain 
numerically compufable. 


The paper is organized as follows. In Secfion|2l we presenf fhe basic fools essentially of complex 
analysis required for fhe sequel. Secfion[3]is devofed to fhe presenfafion of fhe sfrafegy in fhe case of a 
single inclusion, while Secfion|4] deals wifh a finife number of holes separafed by various disfance. 
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2 Basic tools 


In this section, we introduce material which will be used in the following sections. The explicit solution 
to the Laplace problem is well known in the full plane, and we can also find in literature the Green 
function in or outside the unit disk. Conformal mapping is a convenient change of variable for the 
Laplace problem in order to get a formula inside or outside any simply connected compact set. In the 
sequel, we identify and C through (xi, X 2 ) = xi + ix 2 = 5;. 

2.1 Conformal mapping and Green function 

By the Riemann mapping theorem, there exists a unique biholomorphism T* from \ a;® to B{0, 1) 
such that 

T*(oo) = 00 and arg(7~*)^(oo) = 0, 

which reads in the Laurent series decomposition as 

r(z) = /3'z + ^^, VzG^(0,i?r, (2.1) 

fceN ^ 

with /3® G R+ is called the transfinite diameter (or logarithmic capacity) of w®. In the previous equality, 
the radius R is chosen large enough such that tu® C .6(0, R). A consequence of such a decomposition 
is the existence of a constant C > 0 such that: 

W-Tiz) - < C. (2.2) 

For the Dirichlet problem in an open set Q, the Green function is a function Gq from x to M 
such that 


Gn{x,y) = Gn{y,x), Go(x, y) = 0 if x G cAl, AxGn{x,y) = 6{x - y), 


where S is the Dirac measure centered at the origin. In the full plane, we have G]r 2 (x, y) = ^ In | x — y | 
and outside the unit disk we have 

r ) _ 1 , \x-y\ 

r,_ 


G. 


with the notation 


y = 


2vr |x-y*||y|’ 

y 

\yr 


Thanks to the conformal mapping T* we deduce the Green function in the exterior of tu* 


1 


G^^.-{x,y) = — In 


\T\x)-T\y)\ 


(2.3) 


2 tt |T®(x)-T*(y)*||T®(y)r 

Inside Dq, let us introduce the function Q defined by 

Q[^p\{x) := I 5„Gn(,(x,y)(/7(y) dcr(y), Vx G Dq, 

JdQo 

where Gq^ can be defined as in ^^ 2 ^^ replacing T* by a biholomorphism from Dq fo .6(0,1). 
This function safisfies 

(-Ag[(f]=0 in Do, 

[0[(/?]=(/? onODo. 
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Remark 2.1 For any shape uF we can compute numerically T^, T*, only by solving once a 

Laplace problem. This is discussed in y f771 Chapter 16, section 5], We can also make explicit the 
conformal mapping for some geometries. 

When {jF is an ellipse {(x, y) G < c^} with a,b,c> 0, then we find 




and 


V{z) = 


c 


± s/c-'^z^ + - a 


-2 


1-1 


+ b 


-1 


Even if it is not in the geometrical setting of this paper, let us notice that if ui"^ is a segment, then the 
Joukow ski function h{z) = (z + llz)j2 gives an explicit formula for T^. Indeed, h maps the exterior 
of the unit disk to the exterior of the segment [—1,1] x {0}. Then, up to a translation, dilation and 
rotation, we find that T^{z) = z ± s/ z‘^ — 1 (where ± is chosen in order that \T^{z)\ > 1, depending 
on the definition of sends the exterior of the segment to the exterior of the unit disk. 


2.2 Dirichlet problem in an exterior domain 

Even if the main goal of the Green function is to produce an explicit solution of the Laplace problem, 
we only use these functions to prove the following Lemma. 

Lemma 2.2 Let uj be a bounded simply connected open subset o/M^. If F G then the 

boundary value problem 

J—A'I' = 0 

I'll = F on doj, 

admits a unique weak solution 'll in the variational space 

(n-W)M2 + W) 

Furthermore, we have the following properties: 

1. For any n G N, the solution 'll can be decomposed as 

n 

^{X) = Y, ^k{X) + Rn+i (X) (2.4) 

fc =0 

Ofc cos{k9) + sin(fc0) 

where iho ^ constant, 9) = -^- and 

\X\'^\Rn{X)\ < C,,||F||looo^) (2.5) 

for some Cn > 0 independent of F. 

2. Let T be the biholomorphism from onto the exterior of the unit disk (such that T{oo) = oo 

and argT'(oo) = 0), we have 


4-0 = ^o[F] 


d'> = ^/,/(>-)F^n5T(F)d.(r). (2,6) 










3. \I'o = 0 if and only if there exist R> 0 and 'I' E \ Bf), R)) such that 

'I' = Re and 


/ -dz = 0. 
ldB{0,R) ^ 


Proof: The well-posedness in the variational space is a standard result coming from Lax-Milgram 

theorem (see e.g. ifTOll l. 

fdi^\ 

Next, we note that V'I'=(^^jis divergence and curl free, so 

V'I'(z) ;= di^{x,y) - id2'^{x,y) 

is holomorphic because it verifies the Cauchy-Riemann equations. Hence, Vth admits a Laurent series 
decomposition on R(0, R)'^, with R such that uj C B{0,R), and as VT' is square integrable, we deduce 
that 

+ 00 

Ck 


W(z) = ^4’ for all z E R(0, R)" 


k=2 


Obviously, the function 


H-OO 


9{z) :=^ 


k=l 


-Cfc+l 

kz^ 


is a holomorphic primitive of V'h. Decomposing the function g in real and imaginary part 

g{z) = gi{x,y) +ig 2 {x,y), 


we verify that 

= \ {digi + \d2g1) + f ( 5 i 52 + \d292) = digi - id2gi, 
where we have used the Cauchy-Riemann equations on g. Therefore, there exists T'o E M such that 

+ 00 +00 

= + (2.7) 

k=l k=l 

with 

Ok = — , and bk = — 

k k 

This ends the proof of the decomposition of Point [T] We establish (12.51) at the end of this proof. 

We use (12.71 ) to prove the third point. If 'Pq = 0^ th^ri y) = Re g{z) with g E R(C \ B{0, R)) 
and we compute by the Cauchy residue theorem that 


'^{x,y) = ^0 +giix,y) = ^'o + Re 5 ( 2 ;) = T-o + Re 


/ - 

JdB(0,R) ^ 

At the opposite, let us assume that there exist R and 


T' = Re T' and 


= 0 . 


$ E R(C\H(0,R)) such that 


L 


dB(t),R) 


$ , 

- dz = 0. 

z 


5 






We decompose in Laurent series and the condition 


L 


( 96 ( 0 ,it) 


-dz = 0 
z 


reads as do = 0. As 'h = Re 'h, we conclude that 'I'o = Re do = 0- This establishes Point [3] 

Concerning the second statement, we write the expression of the solution 'P in terms of the Green 
function: 

d>{X)= [ dnG{Y,X)F{Y)da{Y), (2.8) 

J duj 

where n is the outgoing normal vector of \ aJ. Thanks to the explicit formula of G, we compute: 




T{Y)-T{X)* \ 
\TiY)-T{X)*\^) 


da(y). 


Indeed, the assumption dio G allows us to apply the Kellogg-Warschawski theorem (see 11211 The¬ 
orem 3.6]) to state that DT~^ is continuous up to the boundary. 

When X —)■ oo, we note that T (X) —)• oo and T(26)* —^ 0, hence it is clear that 


A parametrization of doo can be 


y = 7(d) = r-' 



for d€[ 0 , 27 r]. 


Hence, ( 7 ^(d))'’“ = —DT ^{T{Y))T{Y) where we have used the Cauchy-Riemann equations. We 
compute again by the Cauchy-Riemann equations that 

iy(d)|2 = r{Y)iDr^\T{Y))f DT~\r{Y))r{Y) = det DT~Hr{Y))\r{Y)\^ = 

We deduce that 


n(y) 


iy(d)i 


-^deiDT{Y)DT-^{T{Y))T{Y) = - ^/det DT{Y){DT{Y))-^T{Y) 


= -v'd5:i)f(F)f2^r(r) = 


<£ieLT(r). 

v^detDr(y) 


Hence, we get 


AY ). = -^sMm' = -TdiiimF), 

This allows to conclude: 

^0 = 7^/ F{Y)^detDT{Y) di7(y) 

Jdu 
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An important consequence of this equality is the following estimate: 


|'I'o[-P]| < ||-^||l°°( 9 w)- (2-9) 

Now, we note by the Laurent series decomposition (12.71) that there exists R large enough such that 
||'I'(2f)|| < |'I'o[-^]| + II-^IIl°°( 9 i.j) for all X E B{Q,RY. Therefore, Tf is bounded by 2||F||Loo(g^j) 
outside this ball, and in ^ 8 ( 0 , i?) \ tJ we use the maximum principle to state that 


||'I'(2f)|| < ||'I'||L°°(aa;u9Z3(o,it)) < 2 ||F||loo(9^) VX E 8(0, i?) \a;, 


hence 


l'f'llL°°(R2\tJ) < 2||8||Loo(g^^). 


Next, we consider 8 o > 1 such that a; C 8(0, Rq — 1). Since T' is harmonic, the mean value formula 
implies 

ll^'f'llL°°(aB{o,ito)) < 2||'h||Loo(0(o^R(,_i)c) < 4||8 ||loo(9^). 

Combine this estimate with the Cauchy formulas gives 

1 r 


|Cfe| = 


VY{z)z’^-Uz <48g||F||Loo(a,). 


2 ^ 7 !" JdBiO,Ro) 

Therefore, for any |X| > 28o and any n > 1, we have 

|Cfc+l| _ ^ 

k=n k=n 

In 8(0, 2i?o) \ w, it is clear that 


\xnRn{x)\ < Y. ^ ^ s-Rr'II^-llL-ce^). 


n—1 


|Xn 8 „(X)| < |xnvl/(x)| + ixn^hol + ^ixr-^ < C(n, 8 o)||F||Loo(a<,). 


k=l 


k 


This ends the proof of (I2.5I) . 


Remark 2.3 The previous lemma is still available assuming less regularity on uj than C"^. To write 
a representation formula, it is enough to assume that uj is (7^’" with a finite number of corners with 
openings in (0, 27r) H13\ p. 20]. We also used the maximumprincipe in {B\uJ), this requires that {B\uJ) 
is the interior ball property that is ZJ has the exterior ball property m Chapter 3]. This requires that 
only comers with opening in (0 , tt) could be considered. 

If oj is less regular (with crack, for example), the main difficulty is to establish (12.61) or (12.81) . see /@]/. 
Nevertheless, to justify our construction for any inclusion, we only need (12.41) and (12.51) . In the case 
of domain with cracks, decomposition dH is still valid but it requires more sophistical tools to prove 
estimates like (1231) (see Ml/j cind we do not want to enter in this special feature. 


Remark 2.4 When u is the unit ball B, then we recover the statement of SJ] Lemma 2.1 (3)] and 



^0 = 0 . 


The property of the zero mean value is crucial in /[7]/ in the case of circular inclusions : it implies the 
profile decay at infinity and allows to construct the terms of the expansion. In the present paper, even 
if the boundary condition has the zero mean value, the associated profile does not decay at infinity and 
we have to lift 'ho suitably (see Section^. 
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We end this section by recalling the following classical elliptic estimate (see e.g. 0 Theorem 2.10]). 


Lemma 2.5 Let xq E flo (5 > 0 such that B{xq,6) C flo- Then, for any n E N, there exists 
Cn{S) > 0 such that 




for any harmonic function u (i.e. such that Au = 0 on B{xo, 6)). 

3 One inclusion 

In the case of one inclusion, we omit the index 1 and denote to, x^, %, ■ ■ ■ 

3.1 One iteration 

To deal with Equation (11.31) satisfied by r^, we consider the more general boundary values problem: 


—Avs = 0 in rie, 

< Ve = If on Orio) 

Ve = f OndoJe- 


(3.1) 


At the first order, is approximated by uq = ^ [f ], see (12.31 ). which is the solution of 



Avq = 0 in Qq, 

) = If on dflo- 


In order to gain one order in the remainder of the asymptotic expansion of Us, we need to introduce 
some notations. For any function / E we define F E H^/^(c)t(j) by 


F(X) = fix), VX 



e 


By Lemmathere exists a unique function 'h E solution of 



(3.2) 



(3.3) 



Alp = 0 in \ uje, 

= f on diOe- 


Then we have for any x E \ We 


■= A{x) = ^iX)= F^[F]iX), withX = E 
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The function T' can be decomposed in polar coordinates (see (12.41) ') as 


T'(r, 0) = T'o +0) with "^{r,9) ='^^dk{0)r ^ and G Span(sin(/i;-), cos(A:-)). (3.4) 

A:>1 

Note that 4' = and then satisfies estimate (12.51) . Let us define 

Te{x) := eT= eT{x) , (3.5) 

which maps fhe exferior of to the exterior of the disk 5(0, e). 

Inspired by the case of the ball |T], we define fhe fwo main ingredienfs of our consfrucfion. The 
firsf one is adapfed fo liff a consfanf function on doje- In H/og’ unique solution of fhe problem 

f—Au = 0 inM^\a4, 

1 u = c on duj^ 

is fhe consfanf funcfion v = c. Buf fhis generafes a consfanf ferm on d^Q and fhe unique solution of 

( —AV = 0 in rio, 

11 / = c on Orio 

is fhe consfanf c ifself. Hence, we can’f reduce fhe remainder wifh fhis procedure. 

When u! is fhe unif ball, another way to lift the constant c in \ 04 is to consider the function 
X ^ . For general inclusion lo, we introduce the function 

4(x) := ln|7;(x)|, 

which is a solution of 

f — A 4 = 0 in \ 0 %, 
l 4 = lne on Owg. 

As behaves at infinity like In |x|, this function does not belong to Hj^g, but this allows us to find a 
non frivial harmonic exfension of consfanf info The frace of fhe funcfion is non zero frace on 
fhe oufer boundary OQq and if is described by fhe behaviour af infinify of fhe Riemann mapping (12.11) . 
and one gefs thanks to ( 12 . 21 ) 

114 - ln/3| • -X£|||Loo(gs^g) < Me. (3.6) 

The second ingredient is to correct the main order of the traces on the outer boundary Oflo- For this, 
we consider the function w,^ := ^[In/3| • —x^W (see (12.31) 1 verifying 

{ —Awu, = 0 in rio, 

ruo; = ln/3| • —Xel onOflo. 

Notice that the boundary condition ln/3| • —x^l is on Oflo, thus G ( 7 ^( 120 ) (independently of 
ui). Since Qq is a bounded domain in we have by the maximum principle 

ll'w^c,;||Lo°(no) ^ II • -a^£lllL°°(no) ^ I ln(/3diam(Ho))|, (3.7) 

where diam(c)f 2 o) denotes the diameter of flo- 
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Remark 3.1 The function Wuj encodes two informations: the shape of oj via the transfinite diameter jd 
and the location of the inclusion via x^. 

Unfortunately this function produces a trace on the small inclusion duj^. The main idea of the 
construction is that an appropriate linear combination of the profiles and Wi_j allows to reduce the 
error on both doce and clUo- 


Proposition 3.2 Let be the solution of (13.11) . There is an harmonic function defined on Ug such 

that 


Ve{x) = Vo(x) + J^uj[F] 


X — Xp 


- ^o[i"] + 


Vo{Xs) - ^'o[^] 


■rD£,a;(a:) + ere{x), 


£ J ~‘ ' w^{xf) — Ine 
with vq = Q{<p\ defined in (12.31) . T'o[F] in (12.61) . in (I3.2I) - (I3.3I) and, with Te given in (13.51) . 

:= (-e ~ with ip = In \Te\ and Woj = ^[In /3\ ■ —x^l]. 


(3.8) 


Thus we have 

lk£llL°°(anonaa;e) < C [||T’||l°°(9(.j) + ll¥^llL°°(aoo)] • 

Remark 3.3 Thanks to (13.71) . — In e does not vanish for e small enough. So that (13.81) makes sense. 

Proof: The harmonicity of the remainder follows from the harmonic!ty of vq, F^[F], Wcj and ip. 

In order to establish the L°° estimate, we compute the traces of the function Vp on each componant dujp 
and dQ of the boundary. 

• On the outer boundary dLlo, one has for any x G OLIq 


£rp{x) = <p{x) - 




ip{x) +Foj[F] 


X — Xp 


-T'oT 


vojxp) - T'o 
Wuj{xe) - Ine 


(4(x) - ln/3\x - Xp\) 


-Xp\ , uo(x£) - T'o 


+ 


Wuj{xp) - Ine 


(4(x) - ln/3|x - Xp\) 


with 'h(= Ri) defined in (13.41) and 4^0 = Therefore, one has according fo (13.61) and (12.51) . 

£|k£||L°°(ano) < Ci£\\F\f^(Q^) + ehpM {\vo{xp)\ + |4'o|), wifh 4 = —-—-—. 

'^Luy^e) ins 

Since |'ho| < ll-^llL°°(aa;) (cf- (12.91) ) and |uo(x£)| < |4||L°°(9no) (^Y the maximum principle), we gel 

lk£||L°°(90o) - ^ [II-^IIl°°(9(.j) + (^•£||V^IIl°°(90o)] ' 

• On (he boundary of (he inclusion dujp, one has for x G dojp 

no(x£) - T-o 


£rp{x) = f{x) - 


vq{x) + F 


X — Xp 


-4'o + 


Wuj{xp) - Ine 


(Ine - Wujix)) 


, , .T, , “ ^0 /I ^ 

vo(x) - Vo H-^ 5 —(Ine - Wi^(x)) 

w^(xp) - Ine 


= -(^^o(a:) - ?^o(a;£)) - —r-^('Vc^(a:£) - w^(x)). 

Ins 
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As soon as e is small enough, there exists <5 > 0 such that uj^ C B{x^^ |) C B{xe, 5) C Oq. Using the 
mean value formula and Lemma 1231 we get 

<\x- x^\ I|V?;o|Il-(; 5(3,^,|)) + (I'f^oCa^e)! + |^o|) \x - Xe\ II V^o; 

< C [e||(^||Loo(gf^g) + £^£||-^||L°°(aa;)] ■ 


Remark 3.4 in /[7]/, the coefficients {a, /3) in front of i,. and w^} in (13.81) are uniquely determined 
when we try to reduce the trace boundaries on dflo <^nd doj^- This gives respectively the two equations 
of the following system 

f o; + /3 = 0, 

[alne + I3w^{xe) = -vo{xf) + ^'o, 


so that a = —j3 = 


Vo{xfj - 


■■ The scale hr 


Wioixe) - Ine 
case of circular inclusion in m Relation (2.6)]. 


w (x )-ln£ analogous to those appearing in the 


3.2 Recursive construction 

Let us come back to the initial problems (IL1I) - (IL2I) . As previously, we consider — uq which 

satisfies (11.31) . This problem has the form (13.11) with 


p = 0, f = -Uq. 

Thus the function vq = Q[p\ = 0. Recall that we denote F(X) = f{x) for any x = Xe + eX E doj£. 
Let us consider T'o = ^o[^] determined in Lemma 1231 Applying Proposition 13.21 we have 


u, 


.(x) - uoix) = r^,ix) = F^[F] (^) - - 




and 

Notice that here 


Wcjixe) - Ine 

^ellL°°(9Qonai^e) - C'lko||L°°(aa;A- 


X 0 £^ui{x) + £rl(x), (3.9) 


F^[F]{^) = -F^AMix). 

By applying iteratively Proposition 13.21 with 

f <y9^(x) = Tg (x), Vx E clUo, 

1 F^{X) = Tg (x), Vx = Xe + eX E dUe, 


we build a sequence of functions defined on OUq and F^ defined on du) such fhaf, for any > 1, 


N 

Ue{x) = Uo(x) + 

k=0 


G[p^]{x) + F^[F^] + 


g[/](Xe)-^0[F^] 
Wui{xe) - Ine 




+ e^+^rf+\x), (3.10) 


where 


^£'''^llL°°(anon9tJe) ^ C'IIt 


h^{due) ^ F 


fc+ll 
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Since f G H ^+^(rio) (see ([L2]l), then uq G Hi+t‘((^o) and ||mo||l°°(Qo) ^ c||f||H-i+M(no)- 


An important point to notice is that we obtain estimates of the remainders in the L°°-norm by this 
method. By properties of harmonic functions (namely, by the maximum principle and Lemma [231) . we 
state that the estimates also hold in the energy norm on any compact subset of Qg- Such a restriction 
also appears in the method developed by M. Dalla Riva and R Musolino in IHIH. 

Nevertheless, estimates in the energy norm in the full domain can be obtained following the strategy 
used in 111 m where one decomposes the correctors on homogeneous harmonic functions. Using El 
Proposition 3.2], traces of functions are estimated on the singular boundary diOe. 

Remark 3.5 In ([7]/, the support of^ is assumed to be far away from the inclusions. With this assump¬ 
tion, the first term of the expansion can be simplified since the term In (13.101) 

is of the same order of the remainder e r]. (x) and then can be removed. Indeed, using Lemma 12.51 we 
have 


VxGaUo, -^o[i^°]| 

In the same way, we also have 

I - uo{xe) - < C'£|ko||L°°(no)- 

Therefore we recover the expansion in jUTj: 

/ \ / \ / \ 1/ \ 
Ue{x) = Uo[x) -\ - - -r-^- tVe,cj{x) + 6 r^{x). 


< £||V'Uo||L°°(B(a;o,(5/2)) < Ce\\uo\\L<^(B{xo,5)), 

< e||'Uo||L°°(no)- 


4 N inclusions well separated 


We come back now to the framework of the introduction, that is we consider N inclusions (w*)i<j<Ar 
of size e centered at points x\. We shall consider two cases depending on the limits of the distance 
between centers 4 when e —)> 0 : 

de := min |x* — xl\. 

1. The first case is the fixed limit centers: when e —)• 0, the centers tend to Xq G Qq and there exists 
C > 0 such that 

de>C. 

2. The second case presents a third scale 77 (e) between e and 1 such that 

4 > Cr]{e), 


with a positive constant C. This scale 77 (e) is assumed to satisfy 

r 7 (e) —)• 0 and +00 as e ^ 0. (4.1) 

A typical choice for 77 (e) is e“ with a G [0,1) as made in 121111. 
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Our aim is to apply here the strategy introduced in the case of a single defect. For any 1 < z < iV, 
we associate Wi and tOe^j as we have introduced and in Proposition I3.2l in the case of a single 
inclusior0: 

In 17^1 - ruj with ru* = ^[ln;S*| •-x*|], 

with /3* the transfinite diameter of the conformal map which sends \ a;| onto \ ^8(0, e). We 
aim at constructing an asymptotic expansion of the solution of (ED- 

To identify the second term of the asymptotic expansion, we superpose the contribution of each 
inclusion: 


N 

Ue{x) = Uo{x) + (-8^48] 

i=l 


N 

^0,i) + i{x) + srl{x), 

i=l 


where 


(^) = -J^^iiuoKx) = ^0,* + O (^) 


as X —>■ oo. 


We look for coefficients such that the remainder erl is of smaller order than the first two terms on 
d^le- 

• Let us first consider the boundary dQ^. By construction, we have 

N N 


erl(x) = - Y - ^o,ij - Y = ^ ( ( ^ + ^ 

i=l i=l 

Let 1 < i < A^. We have for any x € dul 

erl{x) = - uoix) - (-uq^x) - T'o,*) - ae,inie,iix) 




j¥=i 
=T'n.,; - 






- E (-^">1^1 (^) - '■'".j) - 


0,2 

Let us notice that rog^j = Ine — ruj on diol- Furthermore, since |x — Xe| = 0{ds) on dcol, we have 

for j / i. 




Thus 


erl{x) =T'o,j - a£,*(lne - Wi{x)) + 0 ( — ) - 


4/-E“«i>" 


sT^ 


— Wj (x 


'For shortness, we replace the index w* by i 
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By a Taylor expansion of the functions Wj, we deduce 


erl{x) - as,i (ine - Wi{xl) + 0{e)) + O 

- + 


- 

=T'o,j - ae,i{ln£ - Wi{xl)) - ^ flej {In/syxl - x^| - Wj{xl)) 
+ ( 1 + max |aej| ] O { — 


d. 


Then, one cancels the leading terms (up to the order e/dg) by solving 


M, 


with 


Me = 


/ ae,l \ 

( ^0,1 \ 


; 1 

\^0,N/ 

xl) 

lnW\xi - 

-Wj{xl) 

Ine — 


^f) 


(4.2) 


In the following, we distinguish several configurations according to the behavior of de and the confi¬ 
gurations of the inclusions. In each case, we prove that the matrix is invertible. Then we deduce 

the existence of coefficients and a good estimate for the remainder r^. 


Remark 4.1 For a single inclusion (N = 1), we recover that = i— i^ee 

in s j 


4.1 First case: N inclusions at distance 0 ( 1 ) 

Let us first assume that de = 0(1): there exists c > 0 such that 

We-xl\ e [c~^,c] , Vf/j. 

Note that ln/3-^|x* — xl\ — Wj{x\) = 0(1), then the matrix reads 

Me = IneXTv + 0(1), 

which is invertible and verifies 

Consequently, the coefficients Oej satisfy 


max I Oe ,• I = O - 

j ' Vine 
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For the asymptotic construction, we do not use the main term of the inverse of otherwise the 
remainder would be in 0(1/In e). Solving exactly (14.21) gives convenient coefficients so that the 
remainder is in 0(e): 


N 


N 


Ue = Uq 


+ + onSn^. 


i=l 


i=l 


In the first sum, there is no interaction but just a superposition of the effect of each inclusion separately. 
The interaction appears in the coefficients but only in the second order term (in 0(1/ In^ e)) as we 
can see by using (14.31 ). We can adapt the construction at any order using the inverse to construct 
a suitable linear combination of the lifting terms and decrease the order of the remainder terms. 


Remark 4.2 Let us consider N = 2. Then is given by 


Ai-^ = 


Ine — W2{xl) 


I — ln/3^|a:g — XjI + uii(a:g) 


— — x‘^\ + W2ixl)'' 

lns-W2ix'^) 


with 


5(e) = (ine — wi{xl)) (ine — W 2 {xl)) — (in /3^\xl — x‘^\ — wi{x‘^)) {\n.0^\x\ — xl\ — W 2 {x\)) . 

If u)^ and uj"^ are unit ball, then j3^ = fP = 1 and we recover the expressions obtained in /[Dp. 211] 
when f is supported far away the inclusions. Indeed, Remark \3 .51 allows to replace T'o.i by —uo{xl) 
and to remove [F] ( ) — T'o.i- 


4.2 Second case: N inclusions at distance 0{r]{e)) 

We assume now that the distance between any two inclusions is of order 77(e): there exists c > 0 such 
that 

r7(e)|x* -x^l G [c"\c] , Vi / j. 

Since ln/3-^|x| — xi\ = lnT 7 (e) + 0(1) for any i j and Wj{x\) = 0(1) for any i,j, the matrix Ai^ 
satisfies 

M, = lneX7V + ln77(e)(4f^- Xtv) + 0(1) = + 0(1), 

where Ai^ = (Ine — lnT 7 (e)) In + lnr 7 (e) TLn and TLn is the square matrix of size N with every 
coefficients equal to 1. Since the rank of at is one, there exists an orthogonal matrix P such that 


Pn = 

A . 

■ 

= p 

0 

o\ 


A • 

• V 


VO 

0 / 


P 


-1 


Thus 


M'i = P 


In e + {N — 1) In 77 (e) 


0 

(Ine - lnr7(e)) In-i 


P 


-1 
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which is clearly invertible for e small. The inverse is 


(_y\ 40 )-l = p ( ln£+(Af-l)liu?(£) 


0 

_ 1 -7" 

Ins—ln?7{e) N—1 


P 


-1 


-In + 


1 


Ine — lnr/(e) N 


1 


Ine + (A^ — 1 ) In 77 (e) Ine —In 77 (e) 


Pn- 


Consequently, is invertible and, by Neumann series, we can remark 


1 


M~ = 

^ Ine I 1 _ 

\ -L 1- 


1 1 
■In + j^ 


In £ 


|_1 + ( A ,_ 1 ) 1 mM 


1 - 


Ine - 


Pn + O 


In^ e 


The asymptotic expansion of is thus given by 


N 

Ueix) = Uoix) + (^Pp[F] - ^'O, 


2=1 


/ ^ 0,1 \ 


+ M 


-1 


/ W£, 1 (X)' 


+ 0 


\^0,N/ \^e,N{x) 


r7(e) Ine 


where we have expressed the coefficients according to the resolution of the system (I4.2I) . 

In order to compute the leading corrector, one has to consider the limit of lneA4“^ when e —0. 
The limit matrix depends on the ratio In 77 ( 5 )/Ine. If this ratio has a finite limit I, then I G [0,1) 
according to (14.11) and 




1 T , 1 

--In + 


Ine \l-l N [1 + {N -1)1 1-1 


Pn ■ 


(4.4) 


In expression (14.41) . we see that the situation is now more complex than in the previous case (14.31) since 
all the defects interact at leading order 0 ( 1 / Ine). 

As an example, if r]{e) = e" and N = 2, then I = a and we have 


Z 

Ueix) = Uoix) + Y {puji[F] - ^ 0 , 


2=1 


+ 


1 -a\ /rD£,i(x) 


(1 - a^) Ine \ V-a 1 / V'I' 0 , 2 / V^e, 2 (a^) 


+ 0 


In^(e) 


If f is supported far away the inclusions, we obtain, using Remark [3751 the expansion of iH Section 
4.2]: 

Mxq) ^^ ^(x)) + O 


Ueix) = Uoix) - 


(1 + a) Ine 


lm(e) 


4.3 More complex geometrical settings 

There are various situations which can be dealt by our approach. In order to show that the method is 
versatile, let us shortly consider two more complex geometrical setting. 
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4.3.1 Two inclusions at distance C>(£“) and the others at distance 0{1) 

We assume N > 3 and, with a possible renumbering, that the first two inelusions are at distanee 0(e") 
whereas the others are at distanee 0(1), namely there exists c > 0 

G [c“^, c] and |x* — G [c“^, c], \/i < j, {i, j) ^ (1,2). 

In this ease, has the expansion 

la 0 

a 1 

0 1 

with inverse satisfying 





We notiee that only the first two defeets internet at the leading order 0(1/ In e) while the interaetion 
involving the other defeets is postponed at order 0(1/ In^ e). 

Remark 4.3 If we use the expansion of we have the weak asymptotic expansion for with a 

worse remainder: 


Ue 


N 

=uq + Y. 

i=l 



^0,1 — a^0,2 ^0,2 — a^o,i \ - 

lne(l —a^) lne(l —a^) ^lne(l —a^) 


tUej + O 



This expansion sheds light the fact that at the first order, there is only an interaction between the 
first two inclusions. 


4.3.2 Particular case with 3 inclusions 


Let us analyze a last situation with three inelusions and three seales. We eonsider 0 < /? < a < 1 and 
c > 0 sueh that 

— Xgl G [c“^, c], — Xgl G [c“^, c], fori = 1,2. 

The matrix Ais is sueh that 


Ale = Ine + 0(1) with 


/l a j3\ 
= [a 1 ^ 

V/? 13 ij 


Computing the determinant of we have 

det Ma,p = (a — l){20^ — a — 1), 
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which does not vanish when 0 < /3 < a < 1 (since the trinomial function a i-)- 2a^ — a — 1 is negative 
on (0,1)). Consequently is invertible and the inverse of the leading term is explicit 


= p 


0+2+1/02+8/32 

0 


V 


0 


0 

2 


a+ 2 —y/ 02+8/32 
0 


0 \ 


1-aJ 


P 


-1 


where P is the orthogonal matrix given by 




V3 1 \ 
-V3 1 

0 - 2 / 


Consequently, there is a complete interaction between the three defects. 
When a = /3, we recover the main term in (14.41) for iV = 3. 
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